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I. INTRODUCTION 



A primary goal of nucleus-nucleus scattering has been to learn about nuclear structure. 
This has become even more critical in recent years, when many groups became very active in 
the investigation of the physics of nuclei far from the stability, mainly using nucleus-nucleus 
scattering processes at intermediate energies (Ei a b ~ 100 MeV/nucleon). The theoretical 
complexity of such collisions has given rise to the use of a number of different approximations. 
The adequate theoretical tool for this purpose is Glauber's multiple-scattering theory jlj|. 
It has long been known both for its simplicity and amazing predictive power. One can find 
copious examples in the literature where the Glauber theory allows for a simple physical 
interpretation of experimental results as well as their quantitative analysis 0,0,0). In fact, 
fragmentation reactions of the type discussed here have already been successfully analyzed 
in the framework of Glauber's theory: in one-nucleon-removal reactions, the momentum 
distribution of the outgoing fragment has been shown to reflect the momentum distribution 
of the nucleon which is removed from the surface of the projectile nucleus However, 
because of complications involving multiple scattering processes in nucleus-nucleus collisions, 
a full Glauber multiple scattering expansion is impracticable. Fortunately, the study of many 
direct nuclear processes, e.g. nucleon knockout, or stripping, elastic breakup (diffraction 
dissociation), etc, are possible using the optical limit of the Glauber theory, in which the 
nuclear ground-state densities and the nucleon-nucleon total cross sections are the main 
input. In fact, this method has become one of the main tools in the study of nuclei far 
from stability |5|]. When departures from the optical limit are observed, multiple nucleon- 
nucleon collisions and in-medium effects of the nucleon-nucleon interaction and nucleon- 
nucleon correlations become relevant. 

Very peripheral collisions, with impact parameters just around the sum of the nuclear 
radii (grazing collisions), or larger, are well established tools for studying nuclear properties 
with intermediate energies and relativistic heavy ion collisions jflBjO]. These collisions lead 
to excitation of giant resonances through both electromagnetic and strong interactions. At 
intermediate energy collisions (Ei a b ~ 100 MeV/nucleon), or higher, the collision time is 
short and the action of the short-range nuclear interaction can excite the surface region of 
the colliding nuclei. This excitation can equilibrate forming a compound nucleus, and/or 
give rise to pre-equilibrium emission or other fast dissipation processes. 



2 



An interesting reaction mechanism in high-energy peripheral nucleus-nucleus collisions 
was sug gested by Feshbach and Zabek 0, Q. This mechanism has been applied in refs. 



was suggested bvr esnt 

H HQ, Q, Q Q 



to the calculation of pion production in heavy ion collisions from 
subthreshold to relativistic energies. It is assumed that pions are produced in peripheral 
processes through the excitation of the projectiles to a A-isobar giant resonance. The results 
of these calculations were compared to inclusive pion production data for incident energies 
from 50 MeV to 2 GeV per nucleon. As emphasized by those authors, this comparison is 
not very meaningful at high energy where peripheral processes are expected to contribute 
very little to the total pion production. However, at subthreshold energies, coherent pion 
production should dominate the cross section. This mechanism is known as the nuclear 
Weizsaecker- Williams method. It works as follows. 

The uncertainty relation associated to the variation of the time-dependent nuclear field 
on a scale Az leads a relation between the energy, AE, and momentum transfer, Ap: 

. _, h hv . h . ^ — i . 

AE~ — = —, Ap-— =^ AE = vAp. 
At Az Az 

The last equation on the right is the dispersion relation of a phonon. For typical situations, 
Az is a few fermis and the nuclear interaction pulse carries several hundred MeV. This 
relation can also be directly obtained from the collision kinematics. Let (Ei, pj) be the 
initial momentum of the projectile and (AE, Ap) the energy-momentum transfer in the 
reaction. One has 

P f = Pi- Ap, E f = Ei- AE. 

From these relations one finds 

P.Ap _ AE _ - (AE) 2 + (Ap) 2 + (Mf - Mj) c 4 



Ei 2Ei 
Neglect the term on the right-hand side, one gets 

AE = v • Ap = vAp z , (1) 

where Ap z is the momentum transfer along the longitudinal direction. 

The above relation can only be satisfied for nuclear excitations of very small momentum 
transfers, even for moderately large energy transfers. This is the case for the excitation of 
giant resonances. Thus, the nuclear interaction in grazing nuclear collisions is an effective 
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tool to probe giant resonances (for a review see, e.g. ref. |17|). For very large impact param- 
eters (larger than the sum of the nuclear density radii) only the electromagnetic interaction 
is present, and eq. Q (with v ~ c) is just the energy-momentum relation of a real photon. In 
act, relativistic Coulomb excitation is another useful tool for investigating giant resonances 
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The phonon-like relation, eq. ^ is also a tool for studying nucleon-nucleon short-range 
correlations. The energy in eq. ^ could hardly be absorbed by a single nucleon since it would 



carry the momentum ~ y / 2mAE, which is appreciably larger than that of eq. ^ However, 
the phonon could be absorbed by a correlated nucleon-pair, which can have large kinetic 
energy and small total momentum, when the nucleons move in approximately opposite 
directions. This mechanism has been exploited by previous authors to study the emission 
of correlated pairs in relativistic heavy ion collisions 25, Remarkably, refs. j^l ^| do 



not treat properly the nuclear absorption at small impact parameters, leading to very large 
cross sections for the emission of correlated pairs in peripheral collisions. 

In many-body physics the word correlation is used to indicate effects beyond mean-field 
theories. In nuclear physics one distinguishes between short- and long-range correlations. 
Nuclear collective phenomena such as vibrations and rotations are known to be ruled by 
lone-range correlations. These effects are relatively well known. Short-r ang e correlations 



is also a subject of intensive studies in nuclear physics (see, e.g. |l8|, lljj, |2C| |2l|, |22j, |23(). 
The sources of short-range correlations are the strong repulsive core of the microscopic 
nucleon-nucleon interaction at short internucleon distances. The nucleon-nucleon interac- 
tion becomes strongly repulsive at short distances. The phase shifts for 1 So and 3 Si are 
positive at low, and become negative at higher energies |30j]. This indicates a repulsive core 
at short distances and attraction at long distances. In the nuclear medium this repulsive 
interaction is strongly influenced by Pauli blocking. The search for nuclear phenomena 
showing short-range correlations effects is one of the most discussed topics in the nuclear 
structure community. For the nuclear reaction community, the importance of Pauli correla- 
tions in high energy nucleus-nucleus collisions has prompted the consideration of effects of 
dynamical short-range correlations. When one treats nucleus-nucleus collisions at high ener- 
gies with an optical phase shift function one can include both the center-of-mass correlations 
and two-body correlations in a straightforward manner to obtain a rapidly converging series 
for the physical observables. 
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It would be proper at this time to look for fingerprints of short-range correlations in high- 
energy collisions involving rare nuclear isotopes. Recent experiments on knockout reactions 
seem to indicate a quenching of the spectroscopic factor relative to shell-model predictions 
in neutron-rich nuclei Q|. This reduction is thought to be a consequence of short-range 
correlations which spread the single particle strength to states with higher energies. In fact, 
systematic studies with the A (e, e'p) reaction have provided ample evidence for this quench- 
ing phenomenon j^J. In this context, two-proton knockout reactions with exotic nuclear 
beams seem to be a promising tool to investigate short-range correlations in neutron(proton)- 
rich nuclei j^J. Indeed, for decades two-proton knockout has been considered a valuable 
tool to study short-range correlations in proton-nucleus and electron-nucleus processes (for 



recent work, see e.g. 



22, 



23j). In high-energy nucleus-nucleus collisions, the phonon mecha- 
nism, proposed by Feshbach and Zabek, is a useful guide for the investigation of short-range 
correlations. 

The plan of this paper is as follows. In this work we treat the effects of short-range 
correlations on heavy-ion scattering at high energies. In Sec. 2 the Glauber formalism for 
diffraction dissociation is reviewed. In section 3 this formalism is shown to lead to the 
same result as the traditional DWBA calculations under the proper conditions. This is an 
important point, as diffraction dissociation and DWBA approaches are commonly referred 
to as distinct reaction mechanisms in the literature. In section 4 the role of absorption and 
Lorentz boosts is discussed. In section 5 the formalism is applied to heavy-ion collisions in 
the presence of two-body correlations, showing the connection with the Feshbach and Zabek 
method. The significance of short-range correlations is further discussed. In sec. 6 the 
formalism is applied to carbon-carbon and n Li + 9 Be collisions. In Sec. 7 some concluding 
remarks are made. 



II. DIFFRACTION DISSOCIATION 

Let us consider high energy scattering, so that the energy transfer in the collision, AE, 
is much smaller than the kinetic energy of the colliding nuclei, E. In most cases, one is 
also interested in processes for which the fragments fly in the forward direction, i.e. we will 
also assume that AO <C 1. In such situations the particle wavefunctions are well described 
by eikonal waves j^, i.e. a plane wave distorted by an interaction, V, so that the S- 
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matrix is given by the simple formula S (b) = exp [— (ijfw) J dZ V {R)] , with v equal to the 
projectile velocity and R = (b, Z) the distance between projectile and target (V is assumed 
to be spherically symmetric). Extending this approach to account for scattering of bound 
particles, the initial and final states are given by 

% = & (r) exp (ik • R) , = <t> f (r) S (b) exp (ik ■ R) , (2) 

where faj (r) are the initial and final probability amplitudes (wavefunctions) that a particle 
in the projectile is at a distance r from the center of mass. The particle's S-matrix, S(b), 
accounts for the distortion due to the interaction. 

For a projectile with two-body structure (e.g. a core+valence particle) 

*i = 4>i (r) exp [i (k c • r c + k,, • r„)] 

V f = <f> f (r)S c (b c ) S v (b v ) exp [i (k' c -r c + k' v - r v )\ , (3) 

where now faj (r) are the initial and final intrinsic wavefunctions of the (core+valence 
particle) as a function of r = i"i — r 2 . The relation between the intrinsic, r, and center 
of mass, R, coordinates is given in terms of the mass ratios (3i = rrii/mp. Explicitly, 

= R+/9 c r and r c = R — (3 v r. The core and valence particle S'-matrices, S c (b c ) and S v (b v ), 
account for the distortion due to the interaction with the target. 

The probability amplitude for diffraction dissociation is the overlap between the two 
wavefunctions above, i.e. 

^4(diff) = J d 3 r c d 3 r v (p* f (r) (r) 5 (z c + z v ) S c (b c ) S v (b v ) exp [i (q c ■ r c + q„ • r v )) , (4) 

where q c = k^, — k c is the momentum transfer to the core particle, and accordingly for the 
valence particle. The above formula yields the probability amplitude that the projectile 
starts the collision in a bound state and ends up as two separated pieces. The S'-matrices, 
S c and S v carry all the information about the dissociation mechanism. The delta-function 
S (Z) in eq. 0] was introduced to account for the fact that the S'-matrices calculated in the 
eikonal approximation only depend on the transverse direction. 

It is instructive to follow another argument to obtain eq. HJ If only the core scatters 
elastically, whereas the valence particle remains in its unaltered plane wave state, the final 
projectile wavefunction is given by 

= <p f ( r ) [1 _ Sc (b c )} exp [i (k' c • r c + k^ • v v )} . (5) 
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The factor [1 — S c (b c )] is the amplitude for elastic scattering of the core. The same 
relation can be applied for the valence particle. The diffraction dissociation occurs 
by subtracting the simultaneous scattering of the core+valence particle, represented by 
[1 — S c (b c )] [1 — S v (b v )] , from the independent scattering of core and the valence particle, 



i.e. 



S(diff) = [1-S C (b c )] [1 - S v (&„)] -[1-S C (b c )] -[1-S V (&„)] = S c (b c ) S v (b v ) - 1. (6) 
The factor (-1) is not relevant because of the orthogonality of the wavefunctions <f>i (r) and 



(r). Using A {dm) 



V^ki.ka 



S(di 



(diflf) 



^k'^k^j with v 9 ki,k 2 equal to plane waves, we 
regain eq. 0] We thus see that diffractive dissociation (or elastic nuclear breakup) arises 
from the momentum transfer to each particle due to elastic scattering, subtracting the 
momentum transfer to their center of mass. 

The cross section for the diffraction process <pi (r) —></>/ (r) is given by 

2 



da = p (E) 



d 3 r c d 3 r v <p) (r) 



(r) 8 (z c + z v ) S c (b c ) S v (b v ) exp [i (q c • r c + q,, • r„)] 



(7) 



where p (E) is the density of final states, p (E) = S (Q z ) d 3 q c d 3 q v / (27r) 5 , where Q = q c + q„ is 
the momentum transfer to the center of mass of the projectile. The delta function accounts 
for the conservation of the longitudinal momentum of the projectile arising from the use of 
eikonal wavefunctions (i.e. no dependence on the longitudinal cm. scattering). 

It is important to notice that the above formula is somewhat different than eq. 8 of ref. 
29! ] . In that reference the coordinates r, R were used from the start. One can transform the 
integral of eq. [7|to those variables. The Jacobian of the transformation is equal to one and 
d 3 r c d 3 r v = d 3 rd 3 R, d 3 q c d 3 q v = d 3 qd 3 Q, where q = /3 c q v — (5 V <\ C is the momentum transfer to 
the intrinsic coordinates of the projectile. Thus, in the coordinates r, R, eq. reduces to 



da 



d 3 qd 2 Q 
(2vr) 5 



d 3 rd 2 b 4>) (r) 0, (r) S c (b c ) S v (b v ) exp [i (q • r + Q • b)] 



(8) 



The above formula reduces to eq. 8 of ref. [29] if one sets ft v — 1 and (3 C = 0. In this equa- 
tion, (pf (r) can be taken as any final state of the projectile. Thus, it is not only appropriate 
to calculate diffraction dissociation, but also diffraction excitation. Diffraction excitation 
occurs when the final state <pf (r) is a bound state. If it is a state in the continuum (diffrac- 
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tion dissociation), then <pf (r) should be set to the unity 1 , since the part of the wavefunction 
given by S C S V exp [i (k' c • r c + ■ r„)] already accounts for the proper wavefunction of the 
projectile. A natural improvement of eq. [7|is to include final state interactions between the 
core and the valence particle in the coordinate dependence of <frf (r) . 

Since I claim here that eq. |S1 can also be used for calculating excitation cross sections, it 
is adequate to show the relation of this approach to the traditional DWBA and semiclassical 
methods for nuclear excitation in nucleus-nucleus collisions. We will see that the latter are 
perturbative expansions of the eq. |H1 



III. DWBA AND SEMICLASSICAL METHODS 



On can factorize the S-matrices defined in section 2 for the interaction of the core and 
valence particle with the target in terms of their phase-shifts 

i 



"OO 

% I 

X 



/OO 
dZV(R) . (9) 
-OO 



hv 

In the weak interaction limit, or perturbation limit, the phase-shifts are very small so that 
S c (b c ) S v (b v ) = exp [i (xc + Xv)} - 1 + % + iXv 

= 1 + TwJ VcT (rc) dZc + Twj VvT dZv ' (10) 

The factor 1 does not contribute to the breakup. Thus, inserting the result above in eq. 
El one obtains 

Apwba) - ^ J d 3 r c d 3 r v 4>* f (r) ^ (r) [V cT (r c ) + V vT (r„)] exp [i (q c ■ r c + q,, • r„)] , (11) 

where the integrals over z c and z v in eq. ITUlwere absorbed back to the integrals over r c and 
after use of the delta-function 5 (z c + z v ). The above equation is nothing more than the 
plane-wave Born-approximation (PWBA) amplitude. However, absorption is not treated 
properly. For small values of r c and r,; the phase-shifts are not small and the approximation 
used in eq. El fails. A better approximation is to assume that for small distances, where 
absorption is important, S c (b c ) S v (b v ) ~ S(b), where the right-hand side is the S'-matrix 



1 Neglecting final state interactions. If final state interactions are important, <fif (r) is the distortion correc- 
tion to the plane wave. 
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for the projectile scattering as a whole on the target. Using the coordinates r and R , and 
defining U int (r, R) = V cT (r c ) + V nT (r n ), one gets 

^(dwba) = hvA<pwBA) - J d 3 rd 3 R (p* f (r) exp [z'q • r] (r) C/ int (r, R)^ (6) exp [z'Q • R] . 

(12) 

In elastic scattering, or excitation of collective modes (e.g. giant resonances), the mo- 
mentum transfer to the intrinsic coordinates can be neglected and the equation above can 
be written as 

T ( dwba) = <X H (R) C (r) \U int (r, R)\ X i+) (R) 4>i (r)> , (13) 

which has the known form of the DWBA T-matrix. The scattering phase space now only 
depends on the center of mass momentum transfer Q. When the center of mass scattering 
waves are represented by eikonal wavef unctions, one has 

X { - ] * (R) (R) ^ S (b) exp [iQ ■ R] . (14) 

This shows that the PWBA and the DWBA are perturbative expansions of the diffraction 
dissociation formula EJ 

In DWBA (or in the eikonal approximation, eq. IT4*|) . b does not have the classical meaning 
of an impact parameter. To obtain the semiclassical limit one goes one step further. By 
using eq. H2]and assuming that R depends on time so that R = (b, z = vt), the semiclassical 
scattering amplitude is given by A { ^) dass) = i J d 2 b a[^v. lasg) (b) exp(z'Q • b), where 



a 



(*-/) (ls\ _ 1 C(ls\ I ^3 



(b) = — S (6) / dtd s r exp (iu if t) <j>} (r) E/ int (r,t)& (r) , (15) 



(semiclass) V 7 ^ 

where eq. H]was used (Q Z Z — Uif t). 

The semiclassical probability for the transition (i —>■ f) is obtained from the above equa- 
tion after squaring it of integrating it over Q. One gets a^~*^ = J d 2 b P^^JHiass) (P)> where 
-^somkLssi ( b ) = a (lraficiass) ( b ) ' with b having now the explicit meaning of an impact pa- 



rameter. Thus, fl(senfidass) (^)> * s ^he semiclassical excitation amplitude. Equation is 
well-known (for example in Coulomb excitation at low energies) except that the factor S (b) 
is usually set to one. In high energy collisions it is crucial to keep this factor, as it accounts 
for refraction and absorption at small impact parameters: ^(fe)) 2 = exp [2x^ imag ^] , where 
^(imag) j g ca i cu i a ted with the imaginary part of the optical potential. The derivation of the 
DWBA and semiclassical limits of eikonal methods can be easily extended to higher-orders 
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R [fm] 



FIG. 1: Schematic diagram showing how the product of the S-matrix and the interaction potential 
(weighted by the wavefunction) would limit the cross section to grazing impact parameters. 

in the perturbation V. The eikonal method includes all terms of the perturbation series in 
the sudden-collision limit. 

IV. ROLE OF ABSORPTION AND OF LORENTZ BOOSTS 

At this point it is interesting to consider the calculation performed by Feshbach and 
Zabek j^. In that work, eq. EE or its equivalent PWBA form, eq. ^2 without a proper 
account of the strong absorption at small impact parameters (described in eq. El by S (6)), 
was used to calculate the total cross section for emission of a correlated nucleon pair in 
peripheral collisions with heavy ions. Also, interactions without imaginary parts were used. 
As a consequence, they found extremely large cross sections; ~ 1 barn for 16 0+ 16 collisions 
at energies ~ 1 GeV/nucleon. This is certainly inconsistent with perturbation theory. As 
seen schematically in figure ^ the product of the .S-matrix and the interaction potential 
implies that the reaction occurs in a narrow region at "grazing" impact parameters. The 
width of this region is approximately A ~ 1 — 2 fm. The cross section might be written 
as a ~ 27rA (Rp + Rp) P, where P is the average probability for this reaction to occur 
within the impact parameter interval A, and Rp {Rt) is the projectile (target) radius. For 
light nuclei 27rA (Rp + Rt) — 300 — 600 mb. Thus, the probability P violates unitarity 
(perturbation theory is invalid) if cross sections of the order of 1 b are obtained. 
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Ref. |9| also introduced relativistic corrections to the nuclear potential. This relativistic 
property is most easily seen within a folding potential model for a nucleon-nucleus collision: 




dr' 3 



p T (r') v NN (r - r') 



(16) 



where px (r') is the nuclear density of the target. In the frame of reference of the projectile, 
the density of the target looks contracted and particle number conservation leads to the 
relativistic modification of eq. El so that p? (r') — ► 7p^ (r'_j_, jz'), where rj_ is the transverse 



the relative velocity of projectile and target. But the number of nucleons as seen by the 
target (or projectile) per unit area remains the same. In other words, a change of variables 
z" = 72' in the integral of eq. El seems to restore the same eq. El However, this change 
of variables also modifies the nucleon-nucleon interaction vnn- Thus, relativity introduces 
non-trivial effects in a potential model description of nucleus-nucleus scattering at high 
energies. 

Colloquially speaking, nucleus-nucleus scattering at high energies is not simply an inco- 
herent sequence of nucleon-nucleon collisions. Since the nucleons are confined within a box 
(inside the nucleus), Lorentz contraction induces a collective effect: in the extreme limit 
7 — > 00 all nucleons would interact at once with the projectile. This is often neglected in 
pure geometrical (Glauber model) descriptions of nucleus-nucleus collisions at high energies, 
as it is assumed that the nucleons inside "firetubes" scatter independently. 

Assuming that the nucleon-nucleon interaction is of very short range so that the ap- 
proximation t>7VAr (r — r') = Jo 6 (r — r') can be used, one sees from eq. El that V(r), 
the interaction that a nucleon in the projectile has with the target nucleus, also has sim- 
ilar transformation properties as the density: V (r) — > jV (t±,jz), i.e. V (r) transforms 
as the time-component of a four-vector. In this situation, the Lorentz contraction has no 
effect whatsoever in the diffraction dissociation amplitudes, described in the previous sec- 
tions within the eikonal approximation. This is because a change of variables Z' = jZ in 
the eikonal phases leads to the same result as in the non-relativistic case, as can be easily 
checked from eq. El Of course, the delta-function approximation for the nucleon-nucleon 
interaction means that nucleons will scatter at once, and Lorentz contraction does not in- 
troduce any additional collective effect. This is not the case for realistic interactions with 
finite range. Thus, nuclear structure studied with high-energy nucleus-nucleus collisions is 



component of r' and 7 = (1 — v 2 /c 2 ) 



is the Lorentz contraction factor, with v equal to 
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immensely complicated by retardation effects and is not well understood. 



V. EMISSION OF CORRELATED PAIRS IN PERIPHERAL REACTIONS 

Lets us now consider the emission of correlated pairs in peripheral collisions. The pro- 
jectile is now a three-body system, with notation for the coordinates as shown in figure |21 
Following the same arguments used in section 2, the wavefunction of a three-body projectile 
in the initial and final states is given by 

*i = fa (n, r 2 ) exp [i (K c • r c + ki • n + k 2 • r 2 )] 

= 0/ (n, r 2 ) S c (b c ) S 1 (h) S 2 (62) exp [z (K' c • r c + k' x • r a + k' 2 • r 2 )] , (17) 

where now 0jj(r 1; r 2 ) are the initial and final intrinsic wavefunctions of the correlated 
nucleon-nucleon pair as a function of their intrinsic coordinates r 1; r 2 . Assuming that the 
nucleon mass is much smaller than that of the core, one can replace r c ~ R, where R is the 
center of mass of the projectile. 

Following the same steps as before, a relation similar to eq. |H1 can be obtained for the 
cross section for the energy absorption by a correlated pair (when final state interactions 
are neglected): 

d 3 qid 3 q 2 d 2 Q 
(2vr) 8 

x S (b) S 1 (h) S 2 (63) exp [i ( qi • n + q 2 • r 2 + Q • b)]| 2 , (18) 



da 



J d 3 ri d 3 r 2 d 2 b 0} (n, r 2 ) 0, (n, r 2 ) 



where Q = K' c — K c . If the intrinsic nucleon coordinates are denoted by r[ = T{ — R, one 
has bi = \Jb 2 + rf sin 2 6^ + 2rib sin Oi cos {(j) — (pi). 

The above relation can be used for the emission of the nucleon pair. Neglecting final 
state interactions and assuming that the core is not observed (i.e. integrating over Q), one 
gets 



^1 dj~* ^2 (* (* 
da = 1 6 2 / d 2 b \S (b)\ 2 / d 3 r 1 d 3 r 2 fa (n, r 2 ) S 1 (61) S 2 (b 2 ) exp [i (qi ■ n + q 2 ■ r 2 ) 



( 



9) 



In order to proceed further one needs a model wavefunction for the correlated pair, 



hi (r n , r n /). The wavefunction used will have the form 



(r x , r 2 ) = <p a (n) <pp (r 2 ) / corr (r, r c ) , (20) 
12 



FIG. 2: Coordinates used in text for the three-body projectile interacting with the target, the 
coordinate ri is not shown for simplicity. 

where r = ri — r2, <p a (r) = 4> n ij m (r) are single particle wavefunctions with quantum numbers 
a = nljm, and f corr (r, r c ) is a function for the nucleon pair distance r, which also depends 
on a two-particle correlation parameter r c so that / corr (r, r c ) — > as r c — > 0. The effective 
correlation function / corr (r, r c bthe so-called Jastrow factor [32] , is a statistical average of 
the Pauli correlation function [3l| and the correlation function for the dynamical short-range 
(e.g., hard core) correlation. 

A S argued in r e, 0, the true ^d-state wave function of the nucleus containing 
correlations coincide with the independent particle, or Hartree-Fock wavefunction, for in- 
terparticle distances r > r hea i, where r heal ~ 1 fm is the so-called "healing distance". This 
behavior is a consequence of the constraints imposed by the Pauli principle. Nucleons are 
kept apart at short-distances, while for distances beyond several Kp 's there is little effect. 
Consequently, nucleon-nucleon collisions at short distances are rare in nuclear matter, and 
because the strongest part of the interaction is at short distances, the effective force be- 
tween the nucleons is much less than in free space. For example, if a nucleon in 16 O felt the 
cumulative sum of 16 nucleon-nucleon potentials, it would feel a potential of ~ 1400 MeV; 
yet empirically it is known that the effective potential felt by the nucleon in the middle of 
the nucleus is only ~ 40-50 MeV deep. 

Although, in general, the correlation function / corr (r,r c ) may depend on the isospin and 
spin quantum numbers of the two-body channel, we will assume for simplicity that it is a 
plain, state independent, Jastrow factor j^. The effects of nucleon-nucleon correlations in 
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nucleus-nucleus collisions have also been studied in several works. For example, in ref. 
short-range correlations were shown to play an important role in nucleon-nucleus collisions 
at intermediate and high energies. 

The two-particle correlation distance, r c , is a combination of four contributions 



r Pauli + r SRD + r PSR + r CM, 



(21) 



where rp au i; is due to Pauli exclusion-principle correlations, tsrd is related to short-range 
dynamical correlations, tpsr is due to a combination of the Pauli and the short-range dy- 
namical term, and tqm is due to center-of-mass correlations (35]. An approximate set of 
expressions for each of these terms is given by 
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(22) 



1 /3 

where A and Kp = (1.5ir 2 p) ~ 1.36 fm -1 are the target number and the Fermi momentum 
of the target nucleus, respectively, b is a short-range dynamical correlation, b ~ 0.4 fm, B 
is the finite-range parameter of the nucleon-nucleon elastic t-matrix, B ~ 0.62 fm 2 (for col- 
lisions around 200 MeV/nucleon), and l c is the effective "correlation length", l c ~ 1.3 A~ 5 ^ 6 
fm. For proton+ 12 C collisions at 200 MeV/nucleon this set of parameters yields, rp au ii ~ 0.3 
fm, Tsrd — 0.01 fm, rpgp ~ 0.0016 fm, Tcm — 0.18 fm, and r c ~ 0.5 fm. This in fact 
overestimates the correlation distance. A more detailed calculation, using the parameters 
B, b and l c from ref. 34] shows that r c has an appreciable dependence on the collision 
energy, as shown in figure El for protons incident on 12 C. Thus, in nuclear reactions, r c can 
vary substantially with the collision energy and with mass numbers. 

The estimates done above show that the main contribution to the correlation distance 
arises from the Pauli principle. Let us assume a correlation function of the form 

/corr (r) = 1 - exp - v \ 9 * . (23) 
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FIG. 3: Dependence of the short-range correlation distance, r c , on the proton energy in the reaction 
p+ 12 C. 

This correlation function implies that the pair wavefunction decreases for small relative dis- 
tances, r = |ri — r 2 | < r c . The correlation function / corr (r, r c ) goes to one for large values 
of r and to zero for r — > 0. In nuclear structure calculations, the effect of correlation, intro- 
duced by the function / CO rr( r , r c ), becomes large when the correlation distance parameter r c 
becomes large, and vice versa. Here, only the effects of short-range correlations are stud- 
ied and it would be manifest in momentum distributions of highly energetic nucleons, as 
discussed in the introduction, and explicitly shown in the next section. It is important to 
notice that the Gaussian correlation function, eq. EH is unrealistic. Indeed, the short-range 
repulsion is at the origin of the decrease of the pair wavefunction for small relative distances. 
At the same time, there will be an increased probability to find the nucleon pair at medium 
internucleon distances. A two-Gaussian parameterization is needed to quantify this well- 
known effect of short-range correlations. For simplicity, only the simple parameterization of 
eq. I2H1 is used in this work. 

Inserting eqs. EH and ESI i n e q- E3 an d integrating over the pair momenta, one gets 



0"SR 



{C 2 S\ {C 2 S) Vjl 
(2j + l) (2j' + l) 2-, 



d 2 b IS 



x Id n d r 2 \(j) n ijm (ri) S 1 (b x ) (j) n > V j> m > (r 2 ) S 2 (b 2 ) / CO rr(r, rJT 



(24) 



The cross section has been averaged over the initial magnetic quantum numbers of the 
nucleons. If the correlation function were equal to the unity, the integrand would be the 
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product of the probabilities to remove an uncorrelated nucleoli, with quantum numbers 
nljm. The later probability is given by J d 3 r \4> n ij m (r) Si {b)\ 2 . 

The spectroscopic factors in eq. |^]have a complex dependence on the angular momenta of 
the nucleon pair. The correlations arising from angular momentum coupling have been stud- 
ied in ref. |37j. Here we will assume a simple combinatorics so that (C 2 ^)^. = n (n — 1) /2, 
where n is the number of nucleons in the valence shell. 

We see from the equations above that the cross section for the emission of a correlated 
pair is smaller than that for the emission of independent particles, since / corr (r, r c ) < 1. 
Most part of the integrand will have / cor r(r, r c ) ~ 1, except for the small region of volume 
Afr 3 , where Af is a number of order of one. Conservative estimates (using r c = 0.3 — 1 fm), 
imply that the cross section for emission of a correlated pair could not exceed 100 mb, in 
contrast to the results obtained in refs. 0, Q] . We will show this for specific reactions in 
the following section. 



VI. RESULTS AND DISCUSSIONS 

The numerical calculations have been carried out for the systems 12 C+ 12 C at 250 
MeV/nucleon and 11 Li+ 9 Be at 287 MeV/nucleon. In both cases, there are some experimen- 
tal data available for two-nucleon removal. This also allows for the study of the influence 
of a halo wavefunction ( 11 Li) in the results. The wavefunctions were calculated by using a 
Woods-Saxon potential with a spin-orbit and Coulomb potential, 

V(r) = U r (r)+U a (r) + U c (r), (25) 

where 

U r (r) = V r (1 + e'T 1 , U s (r) =V,(1- s) (i + e P^ ( 26 ) 

r dr 

Uc (r) is the potential for a uniformly charged sphere with charge Z — 1 (Z, for neutrons) 
and radius Rc, and pi—(r — Ri) ja^ 

For protons in the 1^3/2 orbital of 12 C the separation energy is 15.96 MeV (the two-proton 
separation energy is 27.18 MeV), which can be reproduced with the parameters V r = —57.41 
MeV, V s = -6.0 MeV, R r = R c = R s = 3.011 fm, a r = 0.52 fm and a s = 0.65 fm. 



16 



The reactions and structure of the two-neutron halo nucleus n Li have attracted much 
interest. It is a Borromean system in the sense that although the three-body system, consist- 
ing of 9 Li and two neutrons, forms a bound state, none of the possible two-body subsystems 
have bound states. Hence the stability of u Li is brought about by the interplay of the 
core-neutron and the neutron-neutron interactions, which must lead to a strongly correlated 
wave function with the two neutrons spatially close together. For the calculation here we 
will approximate the n Li ground state by an inert 9 Li core coupled to a neutron pair in a 
(2sl/2) 2 state, although the most probable configuration is an admixture of neutron pairs in 
(2sl/2) 2 , (lpl/2) 2 , and (ld5/2) 2 states 

However, the former assumption allows for 
a simpler calculation of the correlated-pair emission. The potential parameters are adjusted 
to obtain the single-particle wave functions, reproducing the effective neutron separation 
energies. The two-neutron separation energy is 0.3 MeV. From the systematics in Fig. 6 of 
|4ol ]. the estimated 10 Li average excitation energies is 0.2 MeV for the single-particle state. 
Taking this value for two-neutron coupling to the 9 Li core, one arrives at an effective neutron- 
separation energy of 0.5 MeV. This binding energy for the n+ 10 Li system can be reproduced 
with the potential parameters V r = -42.93 MeV,V s = -6.0 MeV, R T = R c = R s = 3.25 fm 
and a r = a s = 0.65 fm. 

The single-particle wavefunctions obtained in this way were used in eq. ESto calculate the 
momentum distributions of the correlated pair. The integrals in eq. 12*41 were performed using 
a method similar to that described in the appendix of ref. I39J]. The S-matrices (and optical 

rn 

potentials) were calculated by using the "t-pp" interaction, as described in refs. 34, 36]. 
This is the same approximation used in ref. [SJ]. 

In heavy ion physics it is common to define a correlation function by means of 

C(qi ' q2)= (d5^)/(^^)' (27) 

where the cross sections in the denominator are for the emission of a single nucleon. An 
accurate measurement of r c requires the measurement of this correlation function for back- 
to-back (or nearly) pair emission. Until now, heavy ion data refer mainly to small relative 
momentum transfers. Data would only be interesting for the present purposes if the trig- 
gering conditions were changed and if special attention was paid to back-to-back emission. 
In the present work, da / dqidq 2 , instead of C (qi, q2), will be used for the study of emission 
of correlated nucleons. 
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FIG. 4: Contour plots of do jdq\dq 2 for the collision 12 C+ 12 C at 250 MeV/nucleon (upper panels) 
and 11 Li+ 9 Be at 287 MeV/nucleon (lower panels), as a function of p\ = hq\ and p 2 = hq 2 ..Th.e 
left (right) panels are for r c = 0.4 (1) fm. The numbers in the plots indicate the cross section, 
da/dq 1 dq 2 , in units of 10~ 5 [mb/(MeV/c)] 2 . 

In figure |U contour plots for da/dq\dq 2 are presented for the collision 12 C+ 12 C at 250 
MeV/nucleon and 11 Li+ 9 Be at 287 MeV/nucleon, as a function of p\ = Hqi and p 2 = hq 2 . 
The nucleons are assumed to be emitted back-to-back, the nucleon 1 at 0° and nucleon 2 at 
180° with respect to the beam axis, respectively. The upper panels are for the C+C collision, 
while the lower panels are for the Li+Be collisions. The left (right) panels are for r c = 0.4 
fm (r c = 1 fm). The numbers in the plot indicate the cross section, da/dq\dq 2 , in units 
of 10~ 5 [mb/(MeV/c)] 2 . One notices a strong correlation between the nucleon momenta, 
resulting from the phonon relationship, eq. ^ . The effect of the phonon dispersion relation 
is to produce a ridge in the cross section. 

To obtain a greater physical insight, I will now use the PWBA approximation as in eq. 
IT21 (with S(b) = 1), so that, instead of the integrals in eq. UHl one needs now to calculate 

T ( pwba) = J d 3 n d 3 r 2 fa (n, r 2 ) [V x {r x ) + V 2 (r 2 )] exp [i ( qi • n + q 2 ■ r 2 )] . (28) 
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Let us assume that the potentials are given by Gaussian functions, i.e. Vip = 
Vi°2 exp(— r 2 2 /A 2 ) and similarly for the wavef unctions, i.e. <p a ^ = N exp (— r 2 2 /A 2 ). 
In this case it is straightforward to perform analytically all the integrals in eq. I2B If one 
further assumes that the correlation distance, r c , is much smaller than the dimensions of the 
uncorrelated wavefunctions and of the potential, i.e. if r c <C A, A, one gets 

3 



-t(PWBA) — I v 1 + V 2 
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x exp 
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where B is the separation energy of the pair and ^4 = (l + A/A)/2. Therefore, the cross 
section is given by 
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(29) 



Now one can easily understand the physics in figure ^by identifying the terms of the above 
equation. The first term is due to conservation of the momentum along the beam direction, 
which yields the dispersion relation, eq. ^ Note that in the derivation of eq. El it was 
assumed that the core recoils with the same momentum, i.e. Qz — oj/v = — (q± z + q^z)- The 
second term is due to elastic scattering of the pair on the target in the direction transverse 
to the beam. In eq. EH1 both the first and the second terms imply that the momentum 
distribution of correlated pair is such that q\ + g 2 = u/v, i.e. q\ ~ — q 2 , as expected for 
small r c . Also according to these terms, the distribution is smeared by the range of the 
independent wavefunctions of the pair, i.e. (<? 2 2 ) — 1/A 2 . However, the last term implies 
a smearing, or spreading, of the momentum distribution by a much larger factor (assuming 
r c <^ A), i.e. (q ,2 2 ) — l/^ 2 - This explains all physics presented in figure EJ The second 
exponential term in eq. I2H1 plays no role in the results presented in figure since it is 
identical to one. 

As discussed above, the location of the ridges in figure |U is a kinematical property of the 
phonon absorption mechanism, which is independent of the collision energy. Thus, it should 
be observable in intermediate energy collisions (Ei a b ~ 100 MeV/nucleon), as well as in 
relativistic collisions. One also observes that the momentum distributions are narrower for 
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FIG. 5: Singles spectra for the momentum distribution of a nucleon due to the correlated-pair 
mechanism. 

correlated-pair emission from a halo nucleus. This is due to the low binding energy, which 
yields an extended wavefunction of the nucleons in the halo. This is also seen from the first 
exponential term of eq. EH1 since the two-proton separation energy for 12 C is 27.16 MeV, 
while the two-neutron separation energy is 0.3 MeV. The effective value of A is much smaller 
for the first case, leading to a larger spreading of the momentum distributions. However, 
the last term in eq. |2H is still the dominant one leading to a small overall effect on the 
momentum distribution, as shown next. 

It would be interesting to try to observe the contribution of the emission of correlated 
pairs in singles spectra. This can be obtained by integrating da/dq\dq2 over one of the two 
nucleon momenta. This is shown in figure El for C+C and Li+Be collisions, using r c = 0.7 
fm. One observes that the peak in the singles spectra occurs at p ~ h/r c , as expected 
from the arguments presented above. This should be visible in the spectra of nucleons from 
knockout reactions as a bump at high nucleon momenta. The position of the bump would be 
a direct reading of the short-range correlation distance. Notice, however, that such a bump 
could not be noticeable because it is superposed to a large background of knockout nucleons 
from stripping reactions. Only by doing a measurement of back-to-back pair emission, this 
signature of short-range pair-correlations could be assessed. 

The total cross section for the emission of correlated pairs arising from short-range corre- 
lations can be calculated from eq. |2U Assuming r c = 0.7 fm, the total cross section for the 
emission of high-energy correlated pairs in C+C collisions at 250 MeV/nucleon is cr corr = 0.61 
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mb. The experimental value for two-proton knockout in this collision is 5.88 ± 9.70 mb. For 
11 Li+ 9 Be at 287 MeV/nucleon the correlated pair cross section is av rr = 4.1 mb. These 
cross sections are much smaller than those obtained in refs. 0, Q]. For reasons which 
were explained in the paragraph preceding eq. EE the results obtained here are much more 
reasonable. These cross sections are also much smaller than those for one-nucleon knock- 
out reactions (see, e.g. ref. P). They are also only one of the contributions (i.e. only 
from diffraction dissociation) of the two-proton removal cross section. Another contribution 
(stripping) has not been considered here. Stripping would not contribute to back-to-back 
nucleon emission, with nearly zero total momentum of the pair, but is responsible for the 
largest part of the two-nucleon knockout total cross section. 

In conclusion, I have shown that when a projectile reacts with a light nuclear target, the 
short-range correlations contribute to the emission of high-energy nucleons which can be 
visible in measurements of back-to-back emission of nucleon pairs. More experiments and 
also the development of a more complete reaction theory are interesting challenges. The 
theoretical results suggest that the pattern and absolute magnitudes of the partial cross 
sections can provide specific information on the detailed nature of the states involved. This 
is particularly important in the case of reactions involving neutron-rich and proton-rich 
nuclei, far from the stability valley, for which only nuclear reactions are presently capable 
of probing their internal structure. The results presented here will be valuable as a guide to 
extend these studies towards drip line nuclei and look for effects which cannot be probed in 
(e, e') scattering due to the lack of experimental facilities of electron scattering on drip line 
nuclei. 

I would like to thank Angela Bonaccorso, Kai Hencken and Ian Thompson for beneficial 
discussions. This work was supported by the U. S. Department of Energy under grant No. 
DE-FG02-04ER41338. 
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